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CN| Abstract. The medium around massive stars is sttongly shaped by the stellar winds. 

T-H Those winds depend on various stellar parameters (effective temperature, luminosity, 

chemical composition, rotation, ...), which are varying as a function of the time. Using 
the wind properties obtained by classical stellar evolution code allows for the multi- 
Mh D hydrodynamical simulation of the circum-stellar medium accounting for the time 

variations of the wind. We present here the preliminary results of the simulation of the 
(— I medium around a fast rotating star. Comparing the results of such simulations with the 

pi I observed properties of the circum-stellar medium will allow to better understand the 

I interactions of the stellar winds with the interstellar medium, and could give hints on 

^ the past history of the mass loss around massive stars, providing constraints for stellar 

models. 

m 

1. Mass loss by fast rotating massive stars 

(N 

^ It is known from a long time that the luminous flux at the surface of a rotating star is not 

^ constant but varies with the colatitude ( von Zeipel 1924). This result was generalised 

to the case of stars presenting a "shellular"-like rotation ^Maeder|1999l ). It follows that 
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the luminous flux is stronger at the poles than in the equatorial regions. According to 
I the so-called CAK-theory (Castor et al. 1975) , the mass loss from hot massive stars is 

^ governed by the interaction of the luminous flux with the spectral lines. One can thus 

J> expect an anisotropic mass loss around such stars (see e.g. 'Maeder' 2002| . 

This anisotropic mass flux around the star modifies the amount of angular mo- 
mentum removed from the surface. However, the account for the anisotropic winds 
in stellar evolution simulations produces only minor diff'erences compared to models 
where this effect is not accounted for ( Georgy et al.|2011 1. The effect of such winds on 



stellar evolution is thus negligible. One can however wonder whether the anisotropic 
stellar winds can modify the circum-stellar medium (CSM), and eventually leave some 
imprints which could be observable around fast rotating massive stars, or during the 
supernova (SN) explosion. 

In order to explore this direction, we study the case of a very fast rotating 20 
star at very low metallicity (Z = 5 • 10~^). As we are interested in the anisotropic winds, 
we want that the model rotates as quick as possible. The initial angular velocity is set 
at Qini = O.VQcrit where Qcnt is the critical angular velocity, defined by the vanishing 



of the equatorial effective gravity under the effect of the centrifugal force ( jMaeder & 
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Figure 1. Left panel: Hertzsprung-Russell diagram of the 20 Mq. Some key- 
points are indicated: the reaching of the critical velocity (Qcrit), the end of the main 
sequence (End H-b), the departure from the critical velocity (end Qdit), the end of 
the central helium burning phase (end He-b) and the end of the central carbon burn- 
ing phase (end C-b). Right panel: Time evolution of the density (black lines) and 
velocity (red lines) of the stellar winds. The solid lines are for the polar regions and 
the dashed ones for the equator. 



Meynet|2000 1. The low metallicity chosen implies a weak mass loss, ensuring a weak 
braking of the stellar surface. Moreover, the stellar model is computed accounting for 
the internal coupling due to magnetic field ( |Maeder & Meynet|20d4 1, also favouring fast 
rotation. The model is computed with the Geneva stellar evolution code ( [Eggenberger 



et al. 2008 ) from the zero age main sequence up to the end of the central silicon burning. 



On Fig. [T] we show on the left panel the Hertzsprung-Russell diagram for this 
model. Due to the low metallicity, the star remains on the blue side during its evolu- 
tion. The critical velocity is reached already during the main sequence. When the star 
becomes redder, its radius inflates, and the surface velocity decreases. The star moves 
away from the critical velocity. The right panel shows the time evolution of the density 
(in red) and of the velocity (in black) of the stellar winds. The velocity of the wind is 
obtained with the relations of Kudritzki & Puis (2000 1. The density is computed with 
the local mass flux (with mass loss rate recipe of Vink et al. 2000, 2001) and the ve- 
locity. The value for the poles (solid lines) are greater than the equatorial one (dashed 
lines), due to the von Zeipel theorem. 



2. Circum-stellar medium simulations 



The simulations of the CSM around the 20 M© model discussed in the previous section 
are performed with the A-Maze hydrodynamical code ( Walder & Folini|2000 1, which 
is a parallelised, adaptive mesh refinement code. The simulations start with a homo- 
geneous medium, with the inner boarding conditions given by the wind data computed 
previously (they are thus time dependant). The chemical composition of the medium 
is also followed (for the main elements, H, He, C, N, and O). One set of simulations is 
performed in full 3D and cover the first few hundred thousands of years. The other set 
is done assuming an axial symmetry, and cover the whole lifetime of the star. 
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Figure 2. Top panel: Density around the star near the middle of the MS. The 
rotation axis is directed along the jc-axis. The physical dimension along the x-axis is 
roughly 30 pc. Medium panel: Same as above panel, but for the fraction of material 
coming from the stellar winds. Bottom panel: Abundance (mass fraction) of '^C 
surrounding the central star. 
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The main results are illustrated on Fig. [2] Let us mention the following points: 

• First of all, the stellar winds carve a strongly anisotropic cavity in the CSM, 
aligned with the rotation axis of the star. The density in the cavity is much smaller 
than in the ISM. Note also that the anisotropic shape of the cavity is not due to 
the 2D treatment, it also develops in the 3D simulations performed. 

• The cavity is progressively filled with the material of the star. After ~ 5.2 Myr, 
the fraction of matter coming from the star in the cavity is roughly 0.4. The 
medium in the cavity is not homogeneous, and is turbulent. 

• Due to a strong mixing throughout the star, the stellar surface is progressively en- 
riched by the products of hydrogen burning: the hydrogen abundance decreases, 
the helium one increases. Due to CNO-cycle, the carbon and oxygen abundances 
decrease while nitrogen abundance increases. These trends are clearly visible in 
the cavity dug by the winds. 

At the end of the nuclear life of a fast rotating massive star, the CSM is thus 
probably strongly shaped by the stellar winds. The explosion of the central star will 
thus occur in a non-isotropic medium. The next step of this work will be to understand 
how such a medium can affect the explosion itself. This is particularly interesting in 
the frame of the long-soft Gamma Ray Bursts progenitors, which are thought to be fast 
rotating, low metallicity stars (see e.g.|Yoon et aH2006|l. 

To improve our simulations, we are performing some additional computations: 
a) To check whether the aspect ratio of the cavity is affected by the Cartesian grid 
used in these preliminary simulations, the A-Maze code was modified, allowing for 
other (AMR) grid geometries, particularly spherical grids, b) Improving the equation 
of state, by adding the effects of partial ionisation of hydrogen and helium due to the 
strong ionising flux of the star, c) Starting from the 2D simulations, we also plan to 
perform 3D simulations of the most interesting phases of the evolution of the CSM, in 
order to check the 2D results and have a more detailed understanding of the CSM. 

References 

Castor, J. I., Abbott, D. C, & Klein, R. I. 1975, ApJ, 195, 157 

Eggenberger, R, Meynet, G., Maeder, A., Hirschi, R., Qiarbonnel, C, Talon, S., & Ekstrom, S. 

2008, Ap&SS, 316, 43 
Georgy, C, Meynet, G., & Maeder, A. 201 1, A&A, 527, A52+. 
Kudritzki, R.-R, & Puis, J. 2000, ARA&A, 38, 613 
Maeder, A. 1999, A&A, 347, 185 

— 2002, A&A, 392, 575 

Maeder, A., & Meynet, G. 2000, A&A, 361, 159. 

— 2004, A&A, 422, 225. 

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2000, A&A, 362, 295. 

— 2001, A&A, 369, 574. 

von Zeipel, H. 1924, MNRAS, 84, 665 

Walder, R., & Folini, D. 2000, in Thermal and Ionization Aspects of Flows from Hot Stars, 
edited by H. Lamers, & A. Sapar, vol. 204 of Astronomical Society of the Pacific Con- 
ference Series, 28 1 

Yoon, S.-C, Langer, N., & Norman, C. 2006, A&A, 460, 199. 



